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Abstract 

We discuss the nonlinear phenomena of irreversible tipping for non-autonomous 
systems where time-varying inputs correspond to a smooth “parameter shift” from one 
asymptotic value to another. We express tipping in terms of properties of local pullback 
attractors and present some results on how nontrivial dynamics for non-autonomous 
systems can be deduced from analysis of the bifurcation diagram for an associated 
autonomous system where parameters are fixed. In particular, we show that there is a 
unique local pullback point attractor associated with each linearly stable equilibrium 
for the past limit. If there is a smooth stable branch of equilibria over the range of 
values of the parameter shift, the pullback attractor will remain close to (track) this 
branch for small enough rates, though larger rates may lead to rate-induced tipping. 
More generally, we show that one can track certain stable paths that go along several 
stable branches by pseudo-orbits of the system, for small enough rates. For these local 
pullback point attractors, we define notions of bifurcation-induced and irreversible rate- 
induced tipping of the non-autonomous system. In one-dimension, we introduce the 
notion of forward basin stability and use this to give a number of sufficient conditions 
for the presence or absence of rate-induced tipping. We apply our results to give criteria 
for irreversible rate-induced tipping in a conceptual climate model example. 


1 Introduction 

In common language, a tipping point or critical transition is a sudden, large and irreversible 
change in output of a complex system in response to a small change of input. There have 
been a range of papers in the applied sciences that use the notion of a tipping point in 
applications such as climate systems PHHUS] or ecosystems [aniiii]. Further work has 
attempted to hnd predictors or early warning signals in terms of changes in noise properties 
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near a tipping point iniEniEsiEiii2H], other recent work on tipping includes |22l 123] • Indeed 
there is a large literature on catastrophe theory (for example [2] and references therein) that 
addresses related questions in the setting of gradient systems. 

There seems to be no universally agreed rigorous mathematical dehnition of what a 
tipping point is, other than it is a type of bifurcation point, though this excludes both 
“noise-induced” [32110] and “rate-induced” [301 [331123] transitions that have been implicated 
in this sort of phenomenon. Other work on critical transitions includes for example im 
that examines the interaction of noise and bifurcation-related tipping in the framework 
of stochastically perturbed fast-slow systems. In an attempt to understand some of these 
issues, a recent paper [3] classihed tipping points into three distinct mechanisms: bifurcation- 
induced (B-tipping), noise-induced (N-tipping) and rate-induced (R-tipping). As noted in [3], 
a tipping point in a real non-autonomous system will typically be a mixture of these effects 
but even in idealized cases it is a challenge to come up with a mathematically rigorous and 
testable dehnition of “tipping point”. 

The aims of this paper are: 

(a) to suggest some dehnitions for B-tipping and R-tipping appropriate for asymptoti¬ 
cally constant “parameter shifts” in terms of pullback attractors for non-autonomous 
systems. 

(b) to show how these dehnitions allow us to exploit the bifurcation diagram of the asso¬ 
ciated autonomous (i.e. quasi-static) system to obtain a lot of information about B- 
and R-tipping of the non-autonomous system associated with a parameter shift. 


(c) to give some illustrative examples of properties that allow us to predict or rule out B- 
and R-tipping in various cases. 

Restricting to simple attractors and asymptotically constant “parameter shifts” allows us to 
obtain a number of rigorous results for a natural class of non-stationary parameter changes 
that may be encountered in natural sciences and engineering. 

The paper is organized as follows: In Section we consider arbitrary dimensional non- 
autonomous systems and introduce a class of parameter shifts 'P(A_,A+) that represent 
smooth (but not necessarily monotonic) changes between asymptotic values A_ and A+. 
For such parameter shifts, we relate stability properties of branches of solutions in the au¬ 
tonomous system with hxed parameters to properties of local pullback point attractors of 
the non-autonomous system with time-varying parameters. This means that the systems we 
consider are asymptotically autonomous [231122] both in forward and backward time. 


After some discussion of (local) pullback attractors, we show in Theorem 2.2 that one 


can associate a unique pullback point attractor with each linearly stable equilibrium X_ for 
the past limit system (i.e. the autonomous system with A = A_). Then, we dehne two types 


of tracking: e-close tracking and end-point tracking. We show in Lemma 2.3 that if there is a 
stable branch of equilibria with no bifurcation points from this X_ to some stable equilibrium 
X_|_ for the future limit system (i.e. the autonomous system with A = A+), then this unique 
pullback point attractor will e-closely track a path of stable equilibria for sufficiently slow 
rates. 
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On the other hand, there may be a number of stable paths connecting X_ to - each 
of these with an image that is a union of stable branches, but possibly passing through 
bifurcation points. In this case, the pullback point attractor can clearly only track at most 
one of these paths. Allowing arbitrarily small perturbations along the path, there will be 
pseudo-orbits that can track any stable path. In Theorem 2A we give such result by stating 
necessary conditions for e-close tracking by pseudo-orbits. 

Section]^ gives some dehnitions for B- and R-tipping that correspond to lack of end-point 
tracking for parameter shifts. In that section we prove several results for one-dimensional 
systems, including the following: 


Certain bifurcation diagrams or parameter ranges will not permit irreversible R-tipping 
(Theorem 3.2, case 1.). 

Other bifurcation diagrams will give irreversible R-tipping for certain parameter shifts 
(Theorem 3.2, cases 2. and 3.). 

There are restrictions on the regions in parameter space where irreversible R-tipping 
can occur (Section |3.3[ ). 


These results use some properties of the paths of attractors. We introduce the notion of 
forward basin stability where at each point in time the basin of attraction of an attractor 
for the associated autonomous system contains all earlier attractors on the path. We prove 
that: 


Forward basin stability guarantees tracking (Theorem 3.2 case 1.). 

Lack of forward basin stability, plus some additional assumptions, gives testable criteria 
for irreversible R-tipping (Theorem 3.2 cases 2. and 3.). 


We also give an example application to a class of global energy-balance climate models in 
Section HI 

In Section we discuss some of the challenges to extend the rigorous results to higher¬ 
dimensional state or parameter spaces, to cases where noise is present and to parameter 
paths that are more general than parameter shifts. 


1.1 The setting: non-autonomous nonlinear systems 

Suppose we have a non-autonomous system of the form 

- = f{x,A{rt)), (1) 

where x G M", A(s) G C^(M, M^), / G C'^(]R"'+‘^, M”), t G M and the “rate” r > 0 is hxed. For 
Xq G M”, f, s G M and starting at x{s) = Xq we write 

:= x{t) 
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to denote the solution cocycle [IB], where x{t) is the corresponding solution of Q. We also 
assume <ht,s(a:o) is dehned for alH > s and xq, and note that it satishes a cocycle equation 

for any t > s > u and xq G M"". Note that, in addition to xq, t and s, depends on the 
shape of the parameter shift A, its rate of change r, as well as on the function / that dehnes 
the system. The key difference from the evolution operator of an autonomous system is that 
a cocycle depends on both the time elapsed {t — s) and the initial time s. 

The model Q can be extended to include stochastic effects, e.g. by considering a stochas¬ 
tic differential equation of the form 

dx = f{x,A{rt))dt + rjdWt (2) 

where Wt represents an n-dimensional Wiener process and r] represents a noise amplitude. 
A random dynamical systems approach [T] gives a solution cocycle that depends also on the 
choice of noise path. 

One can understand a lot about the asymptotic behaviour of Q [or ([^] from the bifur¬ 
cation/attractor diagram of the associated autonomous system 

(3) 

where A is hxed rather than time-dependent. General properties of dynamical systems 
suggest that if r is “small enough” then we can expect, for an appropriate dehnition of 
attractor, that trajectories of Q will closely track a branch of attractors of ([^ if the branch 
continues to be linearly stable over the range of the varying A. In this paper we consider only 
the case of linearly stable equilbrium attractors for ^ where such statements can be made 
rigorous through use of geometric singular perturbation theory ra or uniform hyperbolicity 
[5]. The control theory literature discusses related questions; see for example [HI that allows 
A to depend on state as well. However, close tracking of a continuous path that traverses 
several stable branches of Q cannot always be guaranteed, even if r is small enough. 

If the variation of A along a branch of attractors brings the system to a bifurcation point 
where no stable branches are nearby, we have the ingredients of a bifurcation-induced or 
B-tipping point; no matter how slowly we change the parameter, there will be a sudden 
and irreversible change in the state of the system x on, or near where the parameter passed 
through a dynamic bifurcation mnm- 

However, B-tipping is not the only way to get a sudden irreversible change. Even if 
there is a branch of attractors available for the system to track, a rate-induced tipping or 
R-tipping [3], [3l| can occur; if there is a critical value of the rate r beyond which the system 
cannot track the branch of attractors then it may suddenly move to a different state. This 
happens even though there is no loss of stability in the autonomous system. Nonetheless, the 
bifurcation diagram of the associated autonomous system may give us a lot of information 
and indeed constraints on when R-tipping may or may not happen. Nontrivial behaviour 
that can appear as a result has been recently studied in [Ml EH ES] in terms of canard 
trajectories. 
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Figure 1: Examples of A(s) G V{—2,2). (a) shifts monotonically from —2 to 2 while (b) 
does this shift non-monotonically. 

2 Tracking of attractors for non-autonomous systems 

We consider the behaviour of ([^ for one dimensional (d = 1) “parameter shifts”, i.e. choices 
of A(s) that are C^-smooth (but not necessarily monotonic), bounded and satisfying 

A_ < A(s) < A+ for each (finite) value of s, 

that limit to the ends of the interval, A(s) —?• A± as s —)■ ±oo, where the limiting behaviour 
is asymptotically constant (i.e. dA/ds —)■ 0 as s —>■ ±oo). More precisely. 

Assumption 2.1 We assume 

A(s) G 'P(A_, A+) 

where we consider the subset of (A_, A+)).' 

P(A_,A+) = |a(s) : A_ < A(s) < A+, lim A(s) = A± and lim ^ = 0 

I s^iboo s— >-±00 as 

Figure gives two examples of parameter shifts that we consider, though clearly a parameter 
shift may have an arbitrary number of maxima and minima. 

The set of equilibria for the bifurcation diagram of ([^ is 

X = {{x,\) : /(x,A)=0 and A_ < A < A+} C M” x M. 

Writing df{x, A) to denote the linearization of / with respect to x, we identify the subsets 
of X that are linearly stable, unstable and bifurcation equilibria-. 

A’stab = {( 2 :, A) G A : max{Re(cr(d/(x, A)))} < 0}, 

“Tunst {(x,A) G X : max{Re(cr(d/(x, A)))} > 0}, 

Xhii = {{x, X) E X : max{Re(cr(d/(x, A)))} = 0} 
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where a{M) denotes the spectrum of the matrix M. Note that, in addition to generic 
bifurcations, Ahif may include degenerate bifurcation points and even “bifurcations” that do 
not involve any topological change in the phase portrait, e.g. at A = 0. 

We say X is regular if (i) the set of bifurcation points are isolated (Abif has no accumula¬ 
tion point) and (ii) the only asymptotically stable recurrent sets (attractors) are equilibria. 
The hrst property of a regular X is implied by the Kupka-Smale theorem for flows [T^ Thm 
7.2.13] under generic assumptions, though our assumption is less restrictive and allows the 
possibility of high codimension bifurcations in Abif. The second property of a regular X holds 
for one-dimensional systems, but does not hold in general for higher dimensional systems. 

Assumption 2.2 We assume X is regular. Furthermore, for ease of exposition, we assume 
w.l.o.g. that A_ < A+ throughout and that there are no bifurcations at X±. 

Henceforth, we will denote the closure of U hj U and dehne: 

Definition 2.1 We say a C'^-smooth curve X(A) in X is a branch if it does not intersect 
with Abif except possibly at end points. We say it is a stable branch if it is contained in 
•A^stab ■ 

Definition 2.2 Given a parameter shift A(s) G 'P(A_,A+), we say a continuous curve 
(A(s),A(s)) that limits to some (A±,A±) as s ^ ±oo and whose image lies within X is 
a path. We say it is a stable path if its image lies within Agtab 

Note that branches and paths are defined in terms of solutions to the autonomous system 
(|^ and, to some extent, the shape of the parameter shift A(s), but are independent of the 
rate r. A path can traverse several branches that meet at isolated points in Abif, where the 
path may be nonsmooth. It may also traverse the same branch several times. As illustrated 
in Figure]^ some A(s) can be monotonic functions of s and visit each point on a branch at 
just one value of s, others may have internal minima or maxima and visit the same point 
on a branch at more than one value of s. We mention a minor result that we use in a later 
theorem. 

Lemma 2.1 If X is regular then a stable path traverses a finite number of smooth stable 
branches. 

Proof: Regularity of X in particular means that there are no accumulations of bifurcation 
points Abif. Suppose that a stable path (X(s), A(s)) traverses an inhnite number of smooth 
stable branches and hence passes through an inhnite number of bifurcation points. As the 
path is continuous and bifurcation points are isolated on the path, there must be an ac¬ 
cumulation as s —)■ oo or —oo. However, as X has well-dehned end-points in these limits 
X± = lims_,.±oo -A(s), one of X± is an accumulation point, which gives a contradiction. QED 
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2.1 Local pullback attractors for the non-autonomous system 

We are interested in the influence of parameter shifts A(s) G 'P(A_,A_|_) and rates r on 
pullback attractors of the non-autonomous system Q. Note that the limit r —>■ 0 is singular 
in the sense that r = 0 does not exhibit any change in parameter values with time. In terms 
of the cocycle the following definition is a local version of the pullback point attractor 
discussed in HSl Defn 3.48(ii)] that we have adapted to the context of parameter shifts: 

Definition 2.3 We say a solution x{t) of 0 is a (local) pullback attractor if there is a 
bounded open set U and a time T > 0 such that (i) x(t) G U for all t < —T, and (ii) 

for any y ^ U and t eM. we have 

\^t,t-s{y) - x{t) \ -E 0 as s oo. (4) 

In contrast to m Defn 3.48(ii)] we do not require that the pullback attractor attracts a 
non-autonomous open set U{t). For the asymptotically autonomous case we consider here, 
we can relate each attractor for the past limit autonomous system to a pullback attactor for 
the non-autonomous system via the following theorem, whose proof is given in Appendix A|^ 

Theorem 2.2 For any X_ such that (X_, A_) is in Xstah, there is a unique trajectory Xpb{t) 
of ([ip with Xpb{t) —)■ X_ as t ^ —oo, and this trajectory is a (local) pullback point attractor 
/or Moreover there is a neighbourhood U of X_ and T > 0 such that x{f) is the only 
trajectory remaining in U for all t < —T. 


Theorem 2. 2| shows that any trajectory that limits to a linearly stable equilibrium in the 
past is a pullback point attractor, and these pullback attractors are in one-to-one correspon¬ 
dence with choices of such that (W_, A_) G Tstab- Note that this result is similar to [161 
Thm 3.54(i)], except that their assumption that a: = 0 is invariant means that Xpb{t) = 0. 
Any small enough open neighbourhood U of X_ can be used to define the corresponding 
pullback attractor of Q, in particular by choosing small enough e > 0, the unique pullback 
attractor with past limit X_ can be expressed as 


X 


[A,r,X_] 

pb 


(() := n U 


(5) 


s>0 u>s 


for a given parameter shift A(s) G P(A_,A+), rate r > 0 and initial state X_, where 
denotes a u—dimensional ball of radius e. We do not use ([^ in the following, but note that 
it may be useful as a basis for numerical approximation of the pullback attractor that limits 
to a given X_ for t ——oo. 

^We are indebted to an anonymous referee for suggesting this method of proof. 
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2.2 Tracking of stable branches by pullback attractors 

We define two notions of tracking of branches of stable eqnilibria of (|^ by attractors of Q. 
The strongest notion of tracking we consider is e-close tracking-. 

Definition 2.4 For a given e > 0, we say a (piecewise) solution x{t) of 0 for some r > 0 
and A(s) G V e-close tracks the stable path (X(s), A(s)) if 

\x{t) — X{rt)\ < e (6) 


for all f G M. 

The following is another, somewhat weaker notion of tracking, which we call end-point 
tracking-. 

Definition 2.5 Consider a stable path (X(s),A(s)) from (X_,A_) to (X_|_,A+). We say a 
solution x{t) of 0 end-point tracks this path if it satisfies 

lim x(t) = Xm. 

t^±oo ^ ' 

Note that end-point tracking implies that for all e > 0 there is a T > 0 snch that 0 holds for 
all |t| > T. The following result gives a sufficient condition for close- and end-point tracking 
of a stable branch. 


Lemma 2.3 Suppose that A(s) G 'P(A_, A+) and (X(s), A(s)) is a stable path that traverses 
a stable branch from (X_,A_) to (X+,A+) bounded away from dAif. Then for any e > 0 
there is a vq > t) such that for all 0 < r < vq the pullback attractor x'^pf'''^~\t) satisfies 
e-close tracking of (X(s), A(s)). The pullback attractor also end-point tracks the stable path 
(X(s), A(s)) for all sufficiently small 0 < r. 


Proof: Let (X(s),A(s)) be the stable path whose image is contained within the stable 
branch X(A) bounded away from Ahif and we write X(s) = X(A(s)). Now consider the 
augmented system on given by 


f = /(x.AW) 

ds 

— = r. 
dt 

For r = 0 this has a compact invariant critical manifold 


(7) 


Mo:={(X(A(s)),s) : sGM} 

that is foliated with (neutrally stable) equilibria. Note that the section of Mq with constant s 
limits to X± for s —)■ ±oo. As the branch X(A) is bounded away from any bifurcation points, 
there is a p > 0 such that the normal Jacobian satisfies max{Re(cr((i/(X(s), A(s))))} < —p < 


0 for all s G M. In consequence, the one-dimensional invariant manifold Mq is uniformly 
normally hyperbolic. One can check that Mq is a smooth Riemannian manifold of bounded 
geometry. Applying [I2l Theorem 3.1], there is an ro(e) > 0 such that for any r < tq the 
system Q has a unique normally hyperbolic invariant manifold M^, where the distance 
from Mq to Mr does not exceed e. 

Moreover, this manifold contains a single trajectory parametrized by s G M. In summary, 
for any e > 0 and for all 0 < r < ro(e) there is a trajectory {x{t),rt) G Mj. of (0 with 
\x{t) — X{A{rt))\ < e for all rt G M. 


By Theorem 2.2 there is a neighbourhood U containing X_ such that only one trajectory 
x(t) satishes x(t) G 17 as t —)• —oo. Fixing e > 0 and T > 0 such that x(t) G B^{X{t)) C U 
for all t < —T we can see that the trajectory obtained is the pullback attractor 


x{t) = X. 


[A,r,X:_] 

pb 


it). 


In particular x{t) —)■ X_ as t ^ —oo. 

Now consider the behaviour of x{t) in the limit t ^ oo. If we set x{t) := X_|_ -|- z{t) then 

d z 

— = / (X+ -I- A+) -I- G{z, t) 

where G{z, t) = /(X+ -f- 2 ;, A{rt)) — /(X+ -1- z, A+). Noting that 

d z 

- = df{X^,K)z + G{z,t) + 0{\z\^) 
where the error term is uniform in t. For any e > 0 dehne 


( 8 ) 


:= sup{|G( 2 :, F)| : \z\ < e and t'> t} 


and continuity of / implies that g{t) —>■ 0 as t —)■ cxd. For any e > 0, \x{t) — X(A(rf))| < e/2 
for all r < ro(e/2) and t. The fact that lims^oo 7C(A(s)) = X~^ means there is a T(e) > 0 
such that |X+— X(A(rt))| < e/2 for t > T(e). Hence \z{t)\ < e for all t > T(e). Finally, note 
that from the linear stability of i = df{X^, X+)z, there is an e > 0 such that if \z{T)\ < e 
then z(t) —)■ 0 as t —)■ cx). It follows that x(t) —)■ X+ as t ^ 00 and so the solution x(t) 
end-point tracks the branch. QED 


2.3 Tracking of stable paths by pseudo-orbits 

One might wish to characterize a wider class of stable paths that can potentially be tracked 
by including stable paths that pass through bifurcation points where stable branches meet. 
To do this we have to widen our notion of tracking to consider tracking by pseudo-orbits. 
Let us consider the possible stable paths that may be accessible for a given A. 

^We use k = 2 and a = 0 in the notation of [H] and note that because we restrict to Tjtab there is an 
empty unstable bundle. 
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Definition 2.6 Given a parameter shift A(s) G 'P(A_,A+), we say points (X_,A_) and 
(X+, A+) on Xstah are A-connected if they lie on the same stable path (X(s), A(s)). 


Note that (X_,A_) and (X+, A+) are A(t)-connected if and only if they are A(rf)- 
connected for any r > 0, so the definition does not depend on r. For given (X_,A_) and 
(Af_|_, A+) that are A-connected, it is not necessarily the case that there will be tracking by 
orbits of the non-antonomons system for some r > 0. In particnlar, there is no gnarantee 
that the pnllback attractor from Theorem 2.2 remains close to X{rt), that 


lim \t) = X+ 


t—^OO 


or even that hmt_,.oo is an attractor for the fntnre limit system. Indeed the example 

in Sectionshows that snch a trajectory may not exist: achieving the appropriate branch 
switching may not be possible even on varying r and A. 

As discnssed later in Section 2.4 there may be several stable paths from X_ that end 
at different X+. We show that there will be tracking of any stable path by a pseudo-orbit, 
namely an orbit with occasional arbitrarily small adjustments. Recall the following definition 
of a pseudo-orbit (N.B. there are several possible definitions of pseudo-orbit for flows; e.g. 
see [21]). We say x{t) is an (e, T)-pseudo-orbit [211 of Q on [Ta,Tb] if there is a finite series 
ti G [Ta, Tb] for i = 1,..., M, with ti = Ta and fju = Tb, such that for alH = 1,..., M — 1, 


• ti+l > ti T 

• x{t) is a trajectory of Q on the time-interval 

• \x{ti+i) - limt^t,_^^_ x{t)\ < e. 

We say x(t) defined for t in some (possibly infinite) interval is an (e, T)-pseudo-orbit if it 
is an (e, r)-pseudo-orbit on all finite sub-intervals. For convenience we refer to an (e, 1)- 
pseudo-orbit as an e-pseudo-orbit. 


Theorem 2.4 Suppose that (X_, A_) and (X+, A+) are A-connected by a stable path (X(s), A(s)) 
for some A(s) G 'P(A_, A+). Then for any e > 0 there is a r^ > Q such that for all r < vq 
there is an e-pseudo-orbit x{t) of 0 with 

\x{t) — X{rt)\ < e (9) 


for all t G 


This Theorem (proved below) suggests that if we consider the system perturbed by noise 
of arbitrarily low amplitude, for any noise amplitude there will be an initial condition, a 
rate and a realization of additive noise such that the perturbed system has a trajectory that 
end-point tracks the given stable path. Under additional assumptions, we believe that stable 
paths will be e-close tracked with positive probability. Before proving Theorem 2.4, we give 


a special case whose proof uses similar argument to that of Lemma 2.3 


10 







Lemma 2.5 Suppose that A(s) G 'P(A_,A+) and (X(A(s)), A(s)) traverses a single stable 
branch that is bounded away from A’bif in the interval s G [sa,Sb] (we include the possibility 
Sa = —oo or Sb = oo). Then for any e > 0 there is a ro > 0 such that for all 0 < r < vq 
there is a trajectory x{t) that satisfies e-close tracking of this branch in the sense of for 
t G [sa/r,Sb/r]. 


Proof: This is similar to that of Lemma |2.3| except that instead of s G M we consider 
s G [sa, Sb] and x is not necessarily a pullback attractor for the system, or even unique. Con¬ 
sider the case of Sa and Sb hnite and note that if the closed set {(Af(s), A(s)) : s G [sa,Sb]} 
does not intersect the closed set Xbif then they must be isolated by a neighbourhood. Hence 
there is a p > 0 such that max{Re(cT((i/(X(s), A(s))))} < —p < 0 for all s G S' := [sa,Sb]. 
By considering an augmented system as in the proof of Lemma 2.3 we can apply Fenichel’s 
theorem |T3] to conclude that for small enough r there is a trajectory x{t) that e-tracks the 
branch for rt = s G S. QED 


Proof: (of Theorem 2.4) Let (X(s),A(s)) be the stable path that starts at (X_,A_) and 
ends at (X+, A+). Because X is regular, we can assume that there are at most hnitely many 
bifurcation points on the path (X(s),A(s)). However, it is possible that A(s) may remain 
at bifurcation points for a set of s that may include intervals: there is a hnite set of points 


{si, 




} with 


h < si 


< Sfc < Sfc 


such that A(s) is constant and (X(s), A(s)) G Tbif when s G [s“, sj*"] for some i. Note that the 
generic (transversal) case is s~ = sf but we allow a more general case where the parameter 
can “linger” at bifurcation points for a time. 

We proceed to show that for any e > 0 there is an ro(e) such that it is possible to construct 
an e-pseudo-orbit x{t) that shadows the stable path X{rt) for all r < tq. All discontinuities 
of the constructed pseudo-orbit will be near t G [sj/r, /r]. We show this in detail for 

the hnite interval s G [sh_i,s“]: similar proofs hold for the semi-inhnite intervals (—oo,s]”] 
and [s^,oo) by using uniformity of linear stability near the endpoints X±. Let us dehne 

Pick any e > 0. Continuity of X{s) means we can hnd a 6 such that 


\X{») - Xj\ < i 


( 10 ) 


for all s G [sj — 6, sjj. 

Pick any rj with 0 < rj < 6] Lemma 2.5 implies that there is a rate ro > 0 (depending 
on e and p) and a trajectory x{t) such that for all r < ro there is a trajectory Xj{t) of ([^ 
satisfying 

( 11 ) 


\x{t)-X{rt)\ < - 


on the interval rt G [sA i + p, s ■ — p]. 
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The fact that {Xj,Xj) is a bifurcation point {f{Xj,\j) = df{Xj,Xj) = 0) and f & C‘^ 
means that, by Taylor expansion, for any ,^, (^ > 0 there are M, iV > 0 such that 

|/(x,A)| < M\x - Xj\‘^ + N\X- Xj\ 

whenever \x — Xj\ < ^ and |A — Aj| < C,. Similarly, A(s) G means there is a P > 0 such 
that 

|A(s) - Xj\ < Pd{s, [s7,s+]) 

for all s E [sj — rj, + rj], where d{s, A) = inf{|s — t\ : t E A}. Choosing ( = Prj we have 

\f{x,A{rt))\ < + NPr] (12) 


for all points in 

■= {{x,s) : \x - Xj\ < ^ and s E [sj - T],s^+ r]]}. 
Considering a forward trajectory x{t) of ([^ we have 

\x{t) — x{t')\ < + NPri)\t — t'\ 

for t > T, if {x(t),rt) E Ij{^,ri) between t and t'. Hence by choosing ^ such that 

e 


(13) 

(14) 


e < e < 


and 7] < 


16M ' IbA^P 

(which is possible as long as e < j^) we can ensure that: 

• There is containment C 


From (14), the inequality 


\x{t) — x(u) I < - for all t E [u, u + 2] 


(15) 


holds for any trajectory x(t) such that \x{t) — Xj\ < e for all t E [u,u + 2] where 
{x{t),rt) E 


Because of (15), when rt E [sj —rj, Sj +rj\ we can hnd a hnite number of trajectory segments, 
each of timeTength T > 1, starting at x{u) = X{ru), such that \x(t) — x{u)\ < e/4 for all 
t E [u,u + 2]. 

Hence there is a choice of r; > 0 and rj > 0 such that, if r < rj then there is a pseudo-orbit 
x(t) with the desired property on the interval rt E [sJ — rj, — rj]. More precisely, for any 
rt E S := [sJ — 7], s'j' + rj] we can hnd an increasing sequence of Sj^k and a piecewise solution 
x{t) with = sJ — 7], Sj^M = s'j' + 7] such that Sj^k+i > Sj^k + 1- Combining the estimates 
(lOpT ) and (15) we have for this solution 


\x{t) — X{rt)\ < \x{t) — x{u)\ + \x{u) — X{ru)\ + \X{ru) — Xj\ + \Xj — X{rt)\ 


e e e e 

< tTtTtTt 
4 4 4 4 

= e. 


(16) 
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X 



Figure 2: Construction of a pseudo-orbit x(t) (red) that e-close tracks a stable path X{rt) 
(black) through a bifurcation point. We sketch the case where Sj = -s” = (i.e. the path 

does not linger at the bifurcation point). The short sections near the bifurcation point have 
duration between 1 and 2 time-units so that the width of the smaller intervals of continuity 
in s lie between 1/r and 2/r. If sJ < sh then we choose a pseudo-orbit that remains at Xj 
over this interval. 

for all rf e [sJ — ? 7 , — rf\ except where rt = Sj^k for some k. 

On choosing r < min{rj : j = 1,..., A;} we have a piecewise dehned path x{t) which 
either satishes ([^ at t and ( |I6| holds, or 

I lim x{u) — lim x(u)\ < \ lim x(u) — X(rt)\ < -^ < e. 

By taking the upper limit at any discontinuity, we note that x{t) is an e-pseudo-orbit as 
required. QED 


2.4 An example: dependence on A-connectedness 

Let us consider the system 


— = sin(a:7r)(A(rt) + cos(a:7r)) (17) 

for X, r > 0 and A(s) G V{—2, 2) this is a simple example where we can explore which stable 
paths can be adiabatically tracked for a given A(s). If we note that A_ = —2 and A+ = +2 
then the bifurcation diagram is as in Figure Although there are branches that continue 
through for x = /c G Z, there is an exchange of stability by a combination of pitchfork 
bifurcations. For this system we see that the set of stable branches is connected into one 
component and so for any /c, / G Z there is a stable path from {2k, —2) to (2Z-I-1,2). However, 
if \k — Z| is large then the path will require A(s) with a large number of internal maxima and 
minima. 
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A 

Figure 3: Bifurcation diagram for © ; the solid lines show the stable branches while the 
dashed lines are unstable. Observe that the only possible tracking by a pseudo-orbit that 
starts at (X_,A_) = (0,-2) for a monotonic A(s) as in Figure a) will be the red lines 
shown that hnish at (X+,A+) = (±1,2). However, for A(s) as in Figure [^b) one can also 
track paths by pseudo-orbits that hnish at (X+, A+) = (±3, 2). 


Now consider any trajectory x{t) of ©■ As the lines of equilibria at integer values of x 
remain how-invariant subspaces for any choice of A and r, it follows that 

[a;(0)J < x{t) < [x(0)J ± 1 for all t, 

where [xj is the smallest integer k < x. This means that the pullback attractor from X_ = 0 
cannot track a stable path from (X_,A_) = (0,-2) to (X+,A+) where = ±1,±3,..., 

and A+ = 2 in Figure Indeed, it is easy to see that the pullback attractor in this case is 


simply the constant 
limit system. 


,r,W_] 


= X_ and in all cases this is in fact a repellor for the future 


However, Theorem |2.4| implies that a pseudo-orbit can track such a path - it allows the 
possibility of arbitrarily small adjustments along the way. Starting at (X_,A_) = (0,-2) 
and allowing arbitrarily small adjustments, there are two possible hnishing points for a A(s) 
such as Figure [^a) and four possible hnishing points for a A(s) such as Figure [^b). 


2.5 Bifurcation-equivalent systems 

Since the discussion so far only depends on branches of equilibria, their stabilities and bi¬ 
furcations, many of the properties carry across to systems that have the same bifurcation 
diagram. Suppose that 

? = p( 2 /)/(l/,A(rr)) (18) 
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where y G M"" and p G C'^(]R"',M) is a strictly positive scalar function that is bounded below. 
The associated autonomous system 


^ = p(y)f(yA) 


(19) 


will have the same trajectories in the phase portrait, different magnitude but the same sign 
eigenvalues, and hence the same bifurcation diagram as (|^. For the same reason systems 
Q and (18) will have the same A-connected points X_ and X+. They also have orbits that 
are equivalent in the following (rather weak) sense: 


Theorem 2.6 Consider non-autonomous systems 0 and (T^ - 

(a) Suppose that for some A(rt) G 'P(A_,A+) and r > 0 system 0 has a solution x{t). 
Then there exists a (state-dependent) resealing of time T{t) with T(0) = 0, and a 
yij) = x{T{t)) that solves the equivalent system (18) with A(rr) = A(rT(r)) replaeing 
A(rr). 


(b) Conversely, suppose that for some A(rr) G 'P(A_, A+) and r > 0 system has a 
solution yij). Then, this yfr) defines a (state-dependent) time resealing T{t) with 
T(0) = 0 such that the equivalent system 0 has a solution x{t) = y{T ^(f)) with 
A{rt) = A{rT~^{t)) replacing A{rt). 


Proof: To show (a), dehne a (state-dependent) time rescalling T(r) through a coupled 
system 


dT 

dr 

dy 

dr 


piv), 

p{y) f {y,HrT)), 


with T(0) = 0 and A(rT(r)) = A(rr). Using the assumption on x{t) one can verify that 
{/{t) = x{T{t)) together with a suitable monotonically increasing function T(r) solve the 
coupled system above: 


ft - f (MT(T))A(rT(T))) 

= piy) f (y, A(rT)) = p{y) f{y, A(rr)). 


Hence, |/(r) also solves (18) with A(rr) = A(rT(r)) replacing A(rr). 

To show (6), for a given |/(r) dehne a state-dependent time rescaling through the integral 


T{t) = / g{y{u))du, 


' u=0 
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that satisfies T(0) = 0 and is a monotonically increasing (hence invertible) fnnction of r 
whose first derivative dT/dr = p{y{T)) is bonnded from below. We can now verify that 
x{t = T(r)) = ?/(r) solves ([^ with A(rf) = A (rT“^(f)) replacing A(rf). On the one hand 


dy dx 
dr dr 


dx dT 
dT dr 


pWr))§ 




On the other hand, the assnmption on j/(r) gives 


^ = P{y{r)) /( 2 /(r), A(rr)) = p( 2 /(r)) f{x{t = T(r)), A(rr)) 
= p{y{T)) f {x{t),HrT-\t))) 

= p{y{r)) f . 


Since p{y) > 0, we have that 

f = /(i.A(rt)). 

QED 


In particular, this means that a tracking trajectory in an equivalent system implies that 
there is a tracking trajectory in the first system for a different A. Note that as the transfor¬ 
mation from A to A d oes depend on the trajectory chosen, one needs to be cautious about 
applying Theorem 


2.6 


We also note that varying r can lead to different parameter shifts A. 
Nonetheless, it does allow extension of some results to more general bifurcation equivalent 
systems. 


3 Bifurcation- and rate-induced tipping 

We now propose some definitions and results on bifurcation (B-) and rate (R-) induced 
tipping in the sense of [3]. Recall from [3] that B-tipping is associated with a bifurcation 
point of the corresponding autonomous system, while R-tipping appears only for fast enough 
variation of parameters and is associated with loss of tracking of a stable branch. 


Definition 3.1 Suppose that (X±, A±) G Xstax, and fix A(s) G 'P(A_, A+). Suppose for all 
small enough r > 0 we have 

\imxff^'^-\t) = X+. 


t^OO 


If and X_ are not A-connected then we say there is B-tipping from X_ for this A. 


Definition 3T requires that the system reaches a point on a stable branch that is not 
accessible by following a stable path permitted by the chosen A. The following result shows 
that a bifurcation must indeed be responsible for B-tipping. 
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Corollary 3.1 Suppose there is B-tipping for some A(s) starting at X_. Then there is at 
least one bifurcation point Xb with A_ < Ab < A+ on the stable branch starting at X_. 


Proof: If there is no such bifurcation point, then there must be a stable branch spanning 
from (X_,A_) to (X+,A+). For any parameter shift A(s) G 'P(A_,A+), Lemma 2.3 implies 
that we track this branch for all small enough r > 0 - a contradiction to the assumption that 
there is B-tipping. QED 


We dehne irreversible R-tipping as follows: 

Definition 3.2 Suppose that (X_, A_) G dfstab and fix A(s) G P(A_, A+). We say there is 
irreversible R-tipping from X_ on A if there is a ro > 0 and an X_|_ that is A-connected to 
X_ such that the system end-point tracks a stable path from X_ to X+ for 0 < r < rg but 
not for r = r^, i.e. 


lim X 

t^OO 


[A,r,X_] 

pb 



Y+y^X. 


for 0 < r < ro 
for r = rQ 


Although Theorem 2.4 guarantees e-close tracking by pseudo-orbits for small enough r, an 
increase in r may lead to irreversible R-tipping as above, but then further increase of r may 
lead back to tracking again! Indeed the values of r that give tracking may be interspersed 
by several windows of rates that give tipping; see the example in Section ^ and Figure 
In addition to irreversible R-tipping dehned above, there can be transient R-tipping. 
Recall that our dehnition of end-point tracking makes no assumption that \x{t) — X{t)\ is 
small, except in the limits t —?• ±cx). This means that the system may end-point track the 
stable path for all r > 0, but depart from the path temporarily when r > tq to visit a 
different state/path at intermediate times |3l]- Such transient R-tipping is not discussed 
here. 


3.1 An example with rate-induced tipping for intermediate rates 

Consider a: G M governed by 

dx , ^ . 

— = g{x,A{rt)) (20) 

where A(s) G 'P(A_, A+) and 

g{x, A) = —((a: -|- a -|- hXY + ctanh(A — d)){x — Kj cosh(eA)) (21) 

for a, b, c, d, e and K constants. This example has a particular bifurcation structure as 
illustrated in Figure]^ if we choose 

a = -0.25, b = 1.2, c = -0.4, d = -0.3, e = 3, K = 2. 
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Figure 4: Bifurcation diagrams showing stable (black solid lines) and unstable (black dashed 
lines) equilibria for the system ( 20pl ) with A between A_ = — 1 and A+ = +1. The red lines 
show a numerical approximation of the pullback attractor starting at X_ for parameter shifts 
(22) on varying the rate r. (a) For r < ri the pullback attractor end-point tracks the branch 
through to X_|_. (b) There is a range of rates ri < r < r 2 where the pullback attractor R-tips 
to R+. (c) For r > r 2 the pullback attractor again end-point tracks the branch through to 


Observe that this is a structurally stable bifurcation diagram and so qualitatively robust to 
small changes in any parameter, and indeed to any small changes in g{x,X). For A(s) G 
V{X-, A+) and r > 0 we take solutions of 

^ = -(A-A_)(A-A+). (22) 

Figure illustrates the appearance of critical rates 0 < ri < r 2 such that R-tipping only 
occurs for a hnite range r G (ri,r 2 ). 

3.2 Identifying regions of R-tipping in one-dimensional systems 

In this section we assume one state space dimension {n = 1) throughout and show that 
geometric properties of the bifurcation diagram can be used to infer where R-tipping is or is 
not possible. 

Consider a stable path (X(s),A(s)) for some A(s) G P(A_,A+). For a hxed u G M, 
let B(X(u),A(u)) be the basin of attraction of the stable equilibrium X = X{v) for the 
autonomous system (|^ with parameter A = A(u), i.e. 

B(X(u), A(u)) := {xo : (i(0t(xo, A), X) —)■ 0 as t —)• oo}, 

where 0t(xo, A) is the trajectory of the autonomous system (|^ with initial condition xq. We 
dehne a notion of forward basin stability in terms of A-dependent basin of attraction for the 
attractors of the associated autonomous system on a given stable path. 

Definition 3.3 We say a stable path (X(s),A(s)) for some A(s) G 'P(A_,A+) is forward 
basin-stable if for every u G M the basin of attraction of the stable equilibrium X (v) on the 
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path contains the closure of the set of all earlier positions of this equilibrium along the path: 


{X(m) : M < C B(X(t>), A(i;)) for all w G M. 

Note that forward basin stability: (i) is defined in terms of solutions to the autonomous 
system (|^ and the shape of the parameter shift A(s), but is independent of the rate r, and 
(ii) is a property of a stable path rather than a stable branch traversed by the path. We 
remark that this is a somewhat different notion to basin stability as dehned in [19] : this is a 
quantity that characterises the relative volume of a basin of attraction and is not a property 
of a branch of attractors. 

In the following, we use the notion of forward basin stability to give sufficient conditions 
for irreversible R-tipping in one dimension to be absent (case 1) or present (cases 2, 3). 

Theorem 3.2 Suppose that (X_,A_) and (X+,A+) in A^tab are A-connected by the stable 
path (X(s),A(s)) for some A G P(A_, A+).- 

1. If {X{s), A{s)) is forward basin-stable then x^^’^’^~\t) end-point tracks (X(s), A(s)) for 
all r > 0. Hence in such a case there can be no irreversible R-tipping from X_ on A. 

2. If the stable path (X(s),A(s)) traverses a single stable branch W(A), there is another 
stable brancl^Y {\) with F+ := T(A+) ^ and there are u <v such that 

X(u) gB(F(A(u)),A(u)), 

then (X(s),A(s)) is not forward basin-stable and there is a parameter shift A(s) G 
V{X-, A+) such that there is irreversible R-tipping from X_ on this A. 

3. If the stable path (W(s), A(s)) traverses a single stable branch, and there is a (F+, A+) G 
A^stab with X_ and not A-connected such that 

X_ gB(F+,A+), 

then (X(s),A(s)) is not forward basin-stable and there is irreversible R-tipping from 
X_ for this A. 

Proof: We dehne 


inf X(u), sup Af (u) , 

R>U V^U 

Bfut(n) := 

V>U 


and note that from these dehnitions, Xfut(M) D Xfut(u) and Bfut(M) C Bfut(u) for all u < v. 


^NB this may be defined only on [Aq, A+] where A_ < Aq < A-|.. 


W,t(u) ;= IJ X(u) = 
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For case 1, let us fix any r > 0 and suppose that (X(s),A(s)) is forward basin-stable. 
This implies that X{u) G Xfut(M) nBfut(M) and also that X_ G Bfut(M) for all n G M. For 
any to and xq G Bfut(rto) we dehne x{t) := and claim that 

A(t) := (i(a;(t),Xfut(rto)) 

satishes A(t) —>■ 0 as t —)■ cxd. To see this, note that if A(t) >0 for any t > to then either 

x{t) < inf A(n) or x{t) > sup X{v). 

v>rto v>rtQ 

In either case ;§A(t) < 0 and so A(t) is monotonic decreasing with t. If it converges to some 
Ao > 0 then we can infer there is another equilibrium (not equal to X+) within B(X+, A+) 
for the future limit system: this is a contradiction. Hence for any trajectory x{t) such that 
x{t) —>■ X_ as f —)■ cxD we have 




Since this holds for all to and = X^ we have x{t) —)■ as f —)■ cxd. By 

Theorem 2.2 this trajectory is the pullback attractor and hence the result. 

For case 2 the existence of such u,v implies that X{u) ^ B(X(n),A(n)) and hence 
(A(s), A(s)) is not forward basin stable. We construct a reparametrization 


A(s) := A(a(s)) 


using a smooth monotonic increasing a G M) that increases rapidly from (t(s) = m to 

ct(s) = V but increases slowly otherwise. More precisely, for any M > 0 and > 0 we choose 
a smooth monotonic function cr(s) such that 


cr(s) = S 

for 

s < u, 

1 < £a(s) < M 

for 

U < (j(s) < V 

£a(,) = M 

for 

u + T] < ct(s) < V — r], 

YM = 1 

for 

a{s) > V 


and 


Using the assumption that the stable path traverses a single stable branch, we apply 
Lemma 2.3| to show that for all small enough rates r > 0 we have 

G B(U(n),A(n)). 


Note that by construction. 


X, 


[A,r,X_] 

pb 


(t) = Xpb 


[A,r,X_ 


\t) 


for any t <u/r. Indeed, by picking M large and rj small, we can ensure that 


\x 


[A,r,X- 

pb 


{v/r) -Y{A{v))\ 
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is as small as desired. In particular this means that \v/r) is in the basin of y(A(n)). 

Hence will track Y{A{t)) for t > v/r if r is small enough (recall that s = rt). 

Applying Lemma 2.5 again for the interval s G [n, cxd) means that reducing r, if necessary, 

will give a pullback attractor with —)■ 1+ as t —)• oo. 

For case 3 a similar argument to case 2 implies that the path is not forward basin stable. 
Given any e > 0 there will be a iF > 0 such that if 


I := {X(m') : u' < -K} 


then 


I C[X_- e, A:_ + e] C [X_ - 2e, X_ + 2e] C p| B(y(F), A(n'))- 


v'>K 


By continuity, there is an M such that |/(a:, A)| < M for all x G / and A G [A_,A+]. Now 
pick any 

MK 

r > 2- 

e 

so that for the pullback attractor x{t) := x^^’^’^~\t) we have x{—K/r) G / and 

MK 

\x{K/r) — x{—K/r)\ < 2-< e. 


Therefore 


x{K/r) G Pi B(F(F),A(A)) 


v'>K 


and so x{t) —>■ F+ as f —)■ cxd: there is R-tipping on this path for some 0 < tq < r. 


QED 


3.3 Isolation of B-tipping and R-tipping in one-dimensional sys¬ 
tems 


An interesting consequence of Theorem |3.2| is that rate-induced tipping will appear neither 
for too small a range of parameters, nor on a segment of a branch that lies too close to a 
saddle-node bifurcation. 


Lemma 3.3 If (X_, A_) G Xstah then there is a u > 0 such that for A+ = A_ -|- z/ there is 
a stable branch X(A) on (A_, A+), and no irreversible R-tipping is possible from X_ for any 
AGiP(A_,A+). 


Proof: As X_ is a linearly stable equilibrium, for all small enough <5 > 0 there will be a 
u > 0 such that if A+ = A_ -|- z/ then 

Bs{X_)c IJ B(X(A),A). 

A_<A<A+ 
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Fixing such a 6, continuity of X(A) means that choosing a (possibly smaller) z/ > 0 we can 
ensure 

{X(A) : A_ < A < A+} C Bs{X_) c IJ B(X(A),A). 

A_<A<A+ 

Hence any path on this branch from A_ to A+ will be forward basin stable: by Theorem |3. 2 
case 1, R-tipping is not possible on such a path. QED 


In fact, a stronger statement may be made near a saddle-node bifurcation. Suppose there 
is a saddle-node bifurcation at (Ao,a;o) where the branches exist for A > Aq. The unstable 
and stable branches from the saddle-node bifurcation will be quadratically tangent to the line 
A = Aq in the bifurcation diagram, meaning they will move monotonically in opposite senses 
near to the bifurcation. More precisely, suppose that there is a branch of stable equilibria 
X'^(A) and a branch of unstable equilibria X“(A) < X^(A). Suppose that a F“(A) is the 
nearest branch of unstable equilibria above X^. Then there is a Ai > Aq such that 


max{X“(/r) : Xq < < Ai} < X'^(A) < min{F“(/i) : Xq < fj. < Ai} 


for all Ao < A < Ai. The branch X^(A) is therefore forward basin-stable and by Theorem |3.2 
case 1 there will be no R-tipping for any parameter shifts within this range of parameters: 


this is a stronger conclusion than Lemma 3.3 


As an example see Figure this shows a bifurcation diagram along with three examples 
of choices of (A_,A+). Two of these ranges (R1,R3) will not show R-tipping, while R2 can 
show R-tipping by Theorem 3.2 case 3. 


4 Example: R-tipping in an energy balance climate 
model 


We briefly indicate how we can characterise paths where there will be R-tipping for the 
Fraedrich energy-balance model of global mean temperature: see [3] for more details and an 
example of a path that gives R—tipping. We show how this can be deduced from application 
of Theorem |3.2[ This model for the evolution of global mean surface temperature T is 


^^f{T) = T-'a{-T* + bT^-c), 


(23) 


where the parameters a, b, c, r can be expressed in terms of physical quantities such as the 
solar output, the planetary albedo etc; note that a > 0 and r > 0 from physical arguments. 
We set X = and note that (23) transforms to 


dX dT 2a. ^ . r— 

— = 2T— = -(-X^ + bX - c)VX 
dt dt T 

which (if we are bounded away from X = T = 0) is bifurcation-equivalent to 


(24) 


dX 

dt 


-X^ + bX-c 


(25) 
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Figure 5: Bifurcation diagram for a one-dimensional system with stable branches A and C 
and unstable branch B. For ranges (A_, A+) with A_ < A+ that are illustrated by R1 and R3, 
Theorem 3^ case 1 applies: the paths are forward basin stable meaning there cannot be an 
R-tipping from branch A. For range R2, Theorem 3.2 case 3 implies there can be R-tipping 
from pullback attractors starting on the stable branch A to the stable branch C. 
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in the sense of Theorem |2.6 though we do not need to apply this theorem for the results 


below. Suppose now that 6(A) and c(A) vary smoothly according to some parameter A that 
undergoes a shift (without loss of generality) from 0 to 1. If 6(A)^ > 4c(A) then there are 


two branches of equilibria for (25) 


and hence also for (23); a stable branch 

vTW 


.y-(a) = '!<A + 


4c(A) 


and an unstable branch 


X“(A) = 


6(A) ^biXy - 4c(A) 


These meet at saddle nodes whenever 6(A)^ = 4c(A) and there are no equilibria for 6(A)^ < 
4c(A). The next statement summarizes a number of results that can be proved about this 
system for parameter shifts of A from 0 to 1. 


Theorem 4.1 Consider the system (23) with 6(A) and d{X) depending smoothly on a single 
real parameter A G [0,1]. Suppose that A varies with time by a parameter shift A G V{0, 1). 

(a) If b{X)‘^ > 4c(A) for all A G [0,1] then there is a stable branch from Xf = X®(0) to 
Xf = X*(l) and there is no B-tipping from Xt for any A G ^(0,1). 

(b) Suppose that (a) is satisfied and X^{p) > X'^{i/) for all /r, G [0,1]. Then there is no 
irreversible R-tipping from Xf on any A G ^*(0,1). 


(c) Suppose that (a) is satisfied and there are 0 < /i < z/ < 1 such that X^{p) < X“(z/). 
Then there is a A E 'P{0, 1) for which there is irreversible R-tipping from Xf. 

(d) Suppose that (a) is satisfied and let Xf = X“(0), Xf = X“(l) he the end-points of the 
unstable branch. If Xf < Xf then there is a A E V{0, 1) for which we have irreversible 
R-tipping from Xt_. 


(e) 7/6(0)^ > 4c(0) and there is 0 < p. < 1 such that 6(/i)^ < 4c(p) then there is B-tipping 
for any A E P(0,1) starting at Xf. 


Proof: Conclusion (a) follows using Lemma 2.5 Conclusion (b,c,d) follow as applications of 


Theorem 3.2 cases 1, 2 and 3 respectively. Conclusion (e) follows on verifying Dehnition 3.1 


of B-tipping. 


QED 


Some consequences of this are that (i) if we £x 6 and vary just c, we cannot find an 
R-tipping in this system: this is because in that case 

X“(A) < t < a"(A) 

for all A G [0,1] so that we can apply Theorem |4.1[ b), (ii) similarly we cannot get R-tipping 
if we £x c > 0 and vary just 6. For (i) and (ii), note that it is still possible to get B-tipping 
in either case, (iii) Theorem |4.lK d) can be used to show that the parameter shift considered 
in [3] leads to R-tipping. 
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5 Discussion 


In this paper we introduce a formalism that uses pullback attractors to describe the phe¬ 
nomena of B-tipping and irreversible R-tipping [3] in a class of non-autonomous systems. 
We restrict to a specific type of system, namely those undergoing what we call “parameter 
shifts” between asymptotically autonomous systems and focus on equilibrium attractors. In 
this setting we show in Theorem 2^ that there is a unique pullback point attractor asso¬ 
ciated with each stable equilibrium for the past time asymptotic system, and we can give 
definitions of what it means to “e-close track”, to “end-point track” and to undergo some 
type of “tipping” depending on the details of the system. For one-dimensional systems, we 
obtain a number of results (in particular Theorem 3.2) that give necessary conditions for 
irreversible R-tipping to be present or absent in a system. 

We consider smooth parameter shifts that remain within the end-points of their range and 
that go (possibly non-monotonically) from A_ to A+ for some A_ < A+. It should be possible 
to relax several of these assumptions and give similar results, as long as the parameters have 
well-defined limits for t —)■ ±oo, though this may lead to less intuitive statements. 

We do not consider noise-induced tipping: for this, one needs to go beyond the purely 
topological/geometric setup here, and to discuss probabilities of tipping. We note that 
Theorem |2.4| gives criteria for tracking of stable branches by pseudo-orbits as well as useful 
conditions for the tracking of branches in the presence of low amplitude noise. Indeed, 
similar results can be made about the tracking of unstable branches by pseudo-orbits though 
presumably if these pseudo-orbits are the result of an occasional random process, there will 
not be tracking of unstable branches for long periods of time with any high probability. We 
also do not consider interaction of noise and rate-induced tipping - but see for example [28] . 

Our results rely on an assumption of simple dynamics in state space - in particular that 
the limiting behaviour of the asymptotically autonomous systems are just equilibria, and (for 
the later results) that the state space is one-dimensional and so well-ordered. Although many 
of the phenomena discussed will be present in higher dimensional systems (for example [30]). 
such systems will clearly display additional types of nonlinear behaviour. This includes more 
complex branches of attractors, periodic orbits, homoclinic and heteroclinic connections, 
chaotic dynamics and complex topological structures of attractor branches. Nonetheless, 
notions such as forward basin-stability should generalize to give similar results for tipping 
from equilibria in higher dimensions using parameter shifts of the type considered here. 

As considered in [3] , one particular type of tipping not present in one-dimensional systems 
may already appear in two dimensional state space where the system has a single globally 
attracting invariant set for all values of the parameters. Such a case of transient R-tipping, 
associated with a novel type of excitability threshold in [HU EH ES], cannot be explained 
simply in terms of the branches of the bifurcation diagram even if we include all recurrent 
invariant sets. Instead it is due to the richer geometry of attraction that can appear in two 
or more dimensional systems. 

We end with a brief discussion of an alternative setting that can be used in cases where 
the parameter shift behaviour is determined by an autonomous ODE, for example dX/dt = 
g{X,r) for A G M where g{X_,r) = g{X^,r) = 0 and g{X,r) > 0 for all A G (A_,A+). In this 
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case we can consider an extended system of the form 


d 

dt 



^(A,r) 


As discussed in [26], the presence of irreversible R-tipping can be understood in terms of 
(heteroclinic) connections between saddle equilibria in this extended system and indeed this 
can be used as a tool to numerically locate rates corresponding to R-tipping. By contrast, in 
this paper we make no such assumption and consider the behaviour of the non-autonomous 
system 0 . 
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Appendix A: Proof of Theorem 2.2 

We thank an anonymous referee for suggesting the detailed proof that we reproduce below. 


Proof: (of Theorem 2.2) Firstly, we hx r > 0 and claim there is a unique solution x{t) of 
Q with x{t) —)■ X_ as f —)■ — oo. To see this, let us dehne 


u{e) := sup {\df{x, A{rt)) — df{X_, A{rt))\ : t G M, |a; — X_ | < e} 
6{T) := maxi sup |/(X_,A(rt))|, sup \df{X_,A{rt)) — df{X_,X_)\ 

[_t<-T t<-T 


and note, because f & and /(X_,A_) = 0, that a;(e) —)■ 0 as e —)■ 0 and 5{T) —)■ 0 as 
T —)■ cx). Linear stability means that the eigenvalues of 


A:=df{X_,X_) 

have negative real parts and so there are K > 0 and a > 0 such that \e^^\ < Ke~°‘^ for f > 0 
(in fact, one can choose any a such that max{Re(cr(A))} < —a < 0): see for example [El 
Lemma 3.3.19]). 

Now set X = X_ + y and note that 


y = Ay + h(y, t) 


(26) 


where 

i) ■= fiX- + y, A{rt)) - Ay 

so that 

dh{y, t) =df{X_ + y, A{rt)) - A 

= ldf(X_ + y,A(rt))-df(X_,A(rt)))] + [d/(X_, A(rt))) - d/(X_, A_)] 
+ df(X_,X_) - A 

= ldf(X_ + y,A(rt))-df(X_,A(rt)))] + [d/(X_, A(rt))) - d/(X_, A_)]. 
Hence, for all t < —T we have 


|h(0,f)| <5(T), \dh{yA)\<oji\y\) + 6{T) 


and we seek a solution of (26) such that \y{t)\ < A for t < —T for any T > Tq. We choose 
A > 0 and Tq > 0 such that 

4Ka-^u:{A) < 1, 2Ka-^6{To) < A, and AKa-^6{To) < 1. 

For any continuous y{t) dehned for t < —Tq and \y{t)\ < A we dehne y by 

y{t) = [ e^^~^^'^h{y{u),u)du 
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and note that the reqnired solntion of (26) satishes y = y- Note that 


\m\ < 


+ (a;(A) + 6{T))\y{u)\] du 


' u= — oo 
, -1 


< Ka-^[6{T) + {oj{A) + 6{T))A]. (27) 

and hence ||^|| := sup^^_rp \y{t)\ < A. 

If there are two snch fnnctions 2 / 1,2 (^) one can verify that 

11^1-^211 < - 2 / 2 II 

and so the integral operator is a contraction on the set of continnons y{t) with | 2 /(t)| < A 
for t < —T for the norm || ■ ||. Hence there is a nniqne hxed point snch that \y{t)\ < A for 
t < —T and this is a solntion of (26). Note from (27) with y = y that for t < —T 

|j/(()| < Ka-\to{A) + i(r))||j/|| + Ka-^S(T) < i||j/j| + Kor'HT). 

Hence when T < Tq we have \y{t)\ < 2Ka~^6{T) for t < —T. Hence y{t) —)■ 0 as f —>■ —cxd. 
If we consider U = the A-ball aronnd X_, and for any t < —T dehne 


^pbid') • nu 


(28) 


s>0 u>s 


Then Xphit) is the nniqne trajectory X_ +y{t) for all t < —T and we have \xpbit) — X_| < A 
on f < —T and Xpb —)■ A_ as f —?• — 00 . Note that (28) is well dehned for all t and is invariant 
nnder ([^ - hence it dehnes a nniqne trajectory for all t. 

Moreover, there is a 7 > 0 snch that for any trajectory z{t) with z{t) G U for s < t < —T 
we have 

\z{t) - Xpb{t)\ < -a;pb(s)| 

for s < f < —T. 

To see this, note that \xpb{t) — A_| < A = 2Ka~^5{T) for t < —T, and hence 

\df{xpbit),A{rt)) - A\ < u{A) + S{T). 

It follows from [71 Prop 1] that if Y(t) is the fnndamental matrix solntion of 

y = A{t)y = df{xpb{t),A{rt))y 

then 

\Y{t)Y-\s)\ < Ke-A-^^AX+s{m{t-s) 
ioT s <t < —T. If we choose T large enongh that 

a - K{u:{A) + S{T)) > 0 
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(29) 


and put X = Xpb{t) + y then ([^ can be written 

y = A{t)y + hi{y,t) 

where 


hi{y,t) ■ = 

and so if 11/1 < e then 
where 


f{xpb{t) + y,A{rt)) - f{xpb{t),A{rt)) - df{xpb{t),A{rt))y. 

\hi{y,t)\ < wi(e)ll/| 


a;i(e) := sup{|d/(a; + ?/,A(s)) - ci/(a;, A(s))| : s G M, |a; - X_| < A, |j/| < e}. 


Uniform continuity of df in this region implies that cji —)■ 0 as e —)■ 0. If we £x e small enough 
that 


A:a;i(e) < a - A:(a;(A) + S(T)) 


then it follows from the argument in [HI Thm 9] that if y(t) is a solution of (29) with 
\y{s)\ < e/K for some s < —T then 




for all t with s <t < —T where 7 = a — K{u{A) +6{T)) —Kui{e). Hence \z{t) —Xpb{t)\ —)■ 0 
as t — s —)■ cxD. 


Moreover, note that if y{s) is a solution of (29) where |?/(s)| < e/K for all s < —T we 
have 


for all s, t with s < t < —T. Fixing any t and taking the limit s —)■ —00 implies that y{t) = 0. 
Hence the only solution that satishes \x{t) — Xpb{t)\ < e/K for all t < —T is the pullback 
attractor Xpb- QED 
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